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Abstract. Knowing whether a hadron is formed inside or outside the nuclear 
medium is very important for correctly interpreting jet quenching in heavy-ion 
collisions. The cleanest experimental environment to study the space-time evo- 
lution of hadronization is semi-inclusive DIS on nuclear targets. 2 frameworks 
are presently competing to explain the observed attenuation of hadron pro- 
duction: quark energy loss, with hadron formation outside the nucleus [2,3], 
and nuclear absorption with hadronization starting inside the nucleus [4-6] . I 
demonstrate that the observed approximate A 2 / 3 scaling of experimental data 
cannot conclusively establish the correctness of either energy loss or absorption. 
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In Deep Inelastic Scattering on nuclear targets (nDIS) one observes a suppres- 
sion of hadron production [1] analogous to hadron quenching in heavy-ion colli- 
sions. However, nDIS offers a much cleaner experimental environment to study 
quark fragmentation: the nuclear medium is well known and the multiplicity in the 
final state is low. Knowing how the struck quark propagates in cold nuclear matter 
- most importantly, whether its color is neutralized inside or outside it - is a neces- 
sary prerequisite for correctly using hadron quenching data to extract the unknown 
properties of the hot and dense medium produced at RHIC. Experimental data on 
hadron production in nDIS are usually presented in terms of the "attenuation ratio" 

R%[{z) = N ois dN dz^ I dN dz^ ' i- e -' * ne sm gl e hadron multiplicity on a target 
of mass number A normalized to the multiplicity on a deuterium target. I will only 
analyze the dependence of R\j on the hadron's fractional energy z — E h /v, but 
data also exist binned in the virtual photon energy v or its virtuality Q 2 . 

In the hadron absorption model of Ref. [4] hadronization is assumed to happen in 2 
stages: (i) the struck quark neutralizes its color and forms a so-called "prehadron" 
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h* , and (ii) the observed hadron h is formed. The average formation length of the 
prehadron, (l*)(z, v), and of the hadron, (l h )(z,v), are computed in the framework 
of the standard Lund model. In the HERMES kinematics [1] the prehadron is 
formed well inside the nucleus while the hadron is produced mostly outside. After 
its formation, the prehadron can interact with the surrounding nucleons with a cross 
section er*(^) = 2/3ct' 1 (^) proportional to the experimental hadron-nucleon cross 
section a . The proportionality factor is fitted to n + production data on a Kr target 
at i^boam = 27 GeV 2 [1]. The probability S A h (z, v) that neither the prehadron nor 
the hadron interact can be computed using transport differential equations [4]: 

Si h {z,v) = j d 2 bdyp A {b,y) 

x-y x' x'-x 00 

/P e~7^7 -<x* f dsAp A (b,s) e "W -<*h f dsAp A (b,s) 
dx'Jdx—e * -^-e - (1) 
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where Al = l h — 1* , and p A is the nuclear density. One can recognize exponential 
probability distributions for prehadron and hadron formation. Note also the inte- 
gration over the interaction points (6, y) of the virtual photon 7* with the quark. 
The hadron multiplicity is computed, at leading order in pQCD, as follows: 

1 dN\{z) _ 1 /■ J „ J ..^ , „ 2 ,do* 



N% IS dz 



1 r A n lq 

-X J dxdv J2e}q f (x,Q 2 )—Sf th (z,v)D h f (z,Q 2 ) . (2) 



exp. cuts 



Here a lq and a lA are the lepton-quark and lepton- nucleus cross sections, qf is the /- 
quark distribution function, and D'j its fragmentation function. The model neglects 
hadron elastic scatterings, diffractive hadron production, and resonance production 
and decay. Therefore it is applicable at 0.4 < z < 0.9. 

Energy loss models [2, 3] assume that the struck quark hadronizes well outside 
the medium, and that hadron attenuation is due to medium-induced gluon radi- 
ation off the quark. The quark's energy is reduced from E q = v to E q = v — e, 
where e is the total energy of the radiated gluons, which translates into a modified 
fragmentation function: 

(l-z) 

D){z, Q 2 ; L) = J dAz V(Az ; q, L)_L-D$(_^, Q 2 ) + Po ( q , L) D h f {z, Q 2 ) . 


The "quenching weight" V(Az) is the probability distribution of a fractional energy 
loss Az = e/u, and Pq is the probability of zero energy loss [7]. The quenching 
weight depends on the quark's in- medium path length Lji(b, y) and on the transport 
coefficient q(b,y), defined to take into account non-uniform nuclear density [8]. 
q(0, 0) = 0.5 GeV 2 /fm, is fitted to tt + production on a Kr target analogously to the 
absorption model. Hadron multiplicity is then computed as in Eq. (2), using the 
modified D and S A = 1. Finally, an integration over (b,y) is performed. 



Energy loss vs. hadron absorption 
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Fig. 1. Left: 
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compared to HERMES data at a beam energy of 27 GeV [1]. Right: results of 
the Rm = cA a fit for {He,N,Ne,Kr} at z=0.65 (solid: absorption; dashed: energy 
loss; dotted: data [1]). More z-bins and case studies can be found in [8]. 



In Fig. 1, the 2 models are compared to HERMES data. Both describe well the 
data, and look similar despite the different physics processes. Quite surprisingly, 
this similarity holds up to very heavy targets like Pb. 

A naive argument, often repeated in discussions and seminars on heavy-ion colli- 
sions, is that in first approximation 1 — R M oc A 2 ! 3 in energy loss models because 
the average energy loss (e) oc (L 2 A ) oc A 2 ' 3 . On the other hand, in absorption mod- 
els the survival probability is proportional to the amount of traversed matter, so 
that 1 — R M oc (La) oc A 1 / 3 . Therefore, it is concluded, a simple analysis of the 
A-dependence of R M (or of R\a m heavy-ion collisions) will clearly signal which 
one of the 2 models is correct. 

The above argument is wrong! Where the argument actually fails is for absorp- 
tion models [9]. If the prehadron were produced always at the 7*-quark interaction 
point (i.e., (I*) = 0) then Rm — cA 1 / 3 at all orders in A 1 ! 3 . However, if we al- 
low for a nonzero (I*), its dimensions must be neutralized by the nuclear radius 
Ra, introducing extra powers of A 1 / 3 . Quite generally, if the probability distri- 
bution for the prehadron formation length is finite at zero formation length, then 
R h M oc A 2 ! 3 + O(A), the same power found in energy loss models [4]. This is the 
case for the presented model, as well as for most other absorption models. 

Then, we can hope to distinguish energy loss from hadron absorption by study- 
ing the breaking of the A 2 / 3 law. To this purpose, it was proposed in [4] to select 
a set of targets {A\, A2, . . . , A n }, fix the z bin, and perform a fit of the form 
1 — R M (z) = c(z)A a ( z \ Both c and a must be considered fit parameters for 2 rea- 
sons: first, both contain information on the dynamics of the hadronization process 
[4,8]; second, one can always redefine c in order to absorb a part of a biasing the 
result, so that it is more correct to ask the fit itself what are the correct values of 
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the 2 parameters. The results of the fit are presented in terms of 2cr confidence 
contours in the (c, a) plane, see Fig. 1. This analysis is powerful: it is sensitive 
to model parameters like q [8], and to different physical mechanisms: e.g., partial 
quark deconfinement in nuclei [4]. However, when we apply the fit to the 2 models 
described in the previous section, we have a surprise: energy loss and absorption 
are indistinguishable. The same holds true for all z bins. Increasing the number 
of targets and the span in atomic number does not help in separating the 2 mod- 
els, either, but clearly shows a non negligible breaking of the A 2 ! 3 law at A > 80 [8]. 

In summary, single hadron suppression obeys a A 2 ! 3 law (broken at A > 80) in both 
energy loss and absorption models. Thus, the observed approximate A 2 / 3 scaling 
of experimental data for light nuclei cannot be used as a proof of the energy loss 
mechanism, as is often done. Even the more refined analysis in terms of (c, a) fits 
cannot clearly distinguish the 2 classes of models, though it will help in checking the 
details of the "correct" model after this is established by other means. To answer the 
very important question of whether or not the struck quark starts hadronizing inside 
the nucleus, we need to consider more exclusive observable, like the z-dependence 
of hadron's pp-broadening and Cronin effect [5] or dihadron correlations [10]. 
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